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ABSTRACT 

We study the Photometric Plane (PHP), namely the relation between the 
effective radius r^, the mean surface brightness within that radius <fi>c, and 
the Sersic index n, in optical {R and /) and near-infrared (K) bands for a large 
sample of early-type galaxies (ETGs) in the rich cluster MS 1008 at z = 0.306. 
The PHP relation is logre (1.07±0.06) • logn+(0.219±0.009)- <fi>c + const, 
with an intrinsic dispersion of ^ 32% in ro, and turns out to be independent of 
waveband. This result is consistent with the fact that internal colour gradients 
of ETGs can have only a mild dependence on galaxy luminosity (mass). There 
is no evidence for a significant curvature in the PHP. We show that this can 
be explained if this relation origins from a systematic variation of the specific 
entrop y of ETGs alorig the galaxy sequence, as was suggested from previous 
works i Marauez et alJl2fl01 ). Indeed, considering spherical, non-rotating, one- 
component galaxy models, we find that the specific entropy is exactly a linear 
combination of logrc, </i>c and logn. The intrinsic scatter of the PHP 
is significantly smaller than for other purely photometric relations, such as 
the Kormendy relation and the photometric Fundamental Plane, which is 
constructed by using colours in place of velocity dispersions. The scatter does 
not depend on the waveband and the residuals about the plane do not correlate 
with residuals of the colour-magnitude relation. This implies either that the 
scatter of the PHP does not origin from stellar population parameters or 



* Based on observations collected at European Southern Observatory (ESO ID. 60.A9203, 60.A-9021, 60.O-9025, 66.A-0316). 
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that it is due to a combined effect of such parameters. Finally, we compare 
the coefficients of the PHP at z 0.3 with those of ETGs at z 0, showing 
that the PHP is a valuable tool to constrain the luminosity evolution of ETGs 
with redshift. The slopes of the PHP do not change significantly with redshift, 
while the zero-point is consistent with cosmological dimming of the surface 
brightness in an expanding universe plus the passive fading of galaxy stellar 
populations with a high formation redshift (zf >l~-2). 

Key words: Galaxies: evolution - Galaxies: fundamental parameters - 
Galaxies: clusters: individual: GIG 1008-1224 - Methods: statistical - Tech- 
niques: photometric 



1 INTRODUCTION 



Global properties of early- type galaxies (ETGs), such as luminosities, colours, radii, line 
indices and velocity dispersions, are tightly correlated, implying that stellar population as 
well as dynamical/structural properties of these systems vary smoothly as a function of their 
mass. One of the most well known of such correlations is the Fundamental Plane (FP), which 
is usually expressed as a relation among the effective parameters of ETGs, i.e. the effective 
(half-light) radius and the mean surface brightness wi t hin that radius <u >p, and the 



central velocity dispersion ctq (jPiorgovski &: Davies 



1987 



Dressier et al. 



19871 ). The main 



characteristics of t he FP are its small intrinsic dispersion, in the range of 0.06-0.13 dex 



14 — 3(Wo) in Te (j.Torgensen. Franx fc Kjaergaard 



Pahre. Djorgovski fc de Carvalho 



1998| ). and its tilt, i.e. the deviation of its slopes from those predicted for a virializ ed family 



Tran et al. 



of homologous systems with constant mass-to- light ratios (IBusarello et al.Ml997ri . Several 
works have studied the FP to 2; ~ 1 (e.g. 



2004; 



Wuvts et al. 



2004 and references 



therein). However, due to the high demand of observing time for the measurement of ve- 
locity dispersions, these studies have been based only on small samples of galaxies. For this 
reason, different efforts have been made to construct correlations among purely photometric 
parameters of ETGs, such as the mea n surface brightne ss (or the luminosity) -size relation. 



also known as the Kormendy re . 



(ICaon. Capaccioli fc D'Onofrio 



ation teormendv 



19771 ). and the size-profile shape relation 



1993^ . The main drawback of these relations is that their 



intrinsic dispersion is significantly larger with respect to that of the FP. Since biases on fit- 
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ting coefficients become larger as the dispersion of observed correlations increases, selection 
effects are a crucial issue for the study of purely photometric correlations of ETGs. 



Another interesting correlation among g . 



tometric plane (hereafter PHP, see 



Graham 



obal properties of ETGs is the so-called pho- 



2OO2I GRA02, and references therein), that is 



the correlation among radius, surface brightness and Sersic index (shape parameter) n of 
ETGs. As shown by GRA02, the PHP could have an intrinsic scatter which is comparable 
or slightly larger than that of the FP, therefore making this relation an interesting tool to 
measure galaxy distances and to analyze the properties of galaxies at different redshifts. To 
date, however, a detailed a nalysis of the use of the PHP for studying samples of distant 



galaxies has not been done. 



La Barbera et al 



( 2004 hereafter LMB04) ffistly attempted to 



derive the PHP for ETGs in a cluster at 2 ~ 0.2, finding that a PHP relation seems to exist 
also at this redshift. The other few existing works on the PHP have only analyzed samples of 
galaxies at z ~ 0. Moreover, a straight comparison of results of these works cannot be done, 
since they have derived the PHP in different wavebands, by using str uctural parameters 
defined in d ifferent ways, and for differently selected samples of galaxies 



(2000 



Khosroshahi et al 



2004^ derived the i^-band PHP for ETGs in the Coma cluster and in nearby groups 



using the central surface brightness, Un , the logarithm of the effect ive radius, logrp, and 



the logarithm of the Sersic index, log n. 



Lima Neto et al, 



and 



Marquez et al, 



(20001) 



(hereafter LGM99 and MLCOO respectively) studied the optical PHP of ETGs in clusters 
at z ~ using /xq, the inverse of the Sersic index, z/ = 1/n, and the scale-length of the 
Sersic law (see eq. 7 of LNGOO). All these studies found that the PHP is actually a surface, 
ETGs populating a curved manifold in the space of structural parameters. On the other 
hand, GRA02 derived the PHP in the B band for Virgo and Fornax ETGs, finding that 
ETGs follow a linear relation among log re, </i>e and logn. As noticed by GRA02, since 
the Sersic law fit may not always provide realistic values of the central surface brightness, 
especially for galaxies with high n values, the use of <fi>e should generally be preferred to 
fiQ. We note that this problem becomes even more important at high redshift, where, due to 
seeing effects, the measurement of /iq can require a large extrapolation of the light profile. 

A possible explanation for the existence of a co rrelation among struc tural parameters of 
ETGs has been suggested by LGM99, MLCOO and Marquez et all (|200l[ hereafter MLCOl). 
Using spherical, non-rotating, one-component models of ETGs, these works argued that 
the PHP relation origins from a relation between the mass of ETGs and their specific 



entropy. MLCOO showed that the entropy-mass relation is consistent with that expected in a 
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dissipation-less merging scheme of galaxy evolution, where merging produces a higher level 
of disorder (entropy) in larger galaxies. This interpretation is supported by the fact that the 



PHPs of dwarf a nd normal elliptica. 



normal galaxies (|Khosroshahi et al, 



s app ear to be offset in a direction of larger entropy for 



2004 ). The connection between the PHP and processes 



such as dissipation-less merging makes it even more interesti ng to analyze this relation at 
different redshifts. If ETGs were mostly assembled at 2 > 1 fe.g. lKauffman , one would 

expect that at lower redshift the slopes of the PHP do not change significantly with 2;, while 
its zero-point varies accordingly to the fading of stellar populations. Just as for the FP, this 
would imply that the zero-point of the PHP could be used to measure the formation epoch 
of stellar populations in ETGs and, perhaps, cosmological parameters. Since at different 
redshifts structural parameters are derived in different restframe bands, the evolution with 
z of the PHP can only be addressed on the basis of an accurate knowledge of the wavelength 
dependence of this relation. On the other hand, the dependence of the PHP on waveband 
carries interesting information by itself, depending on how the ratio of structural parameters 
between different wavebands varies along the galaxy sequence. This variation indicates how 
internal colour gradients of ETGs depend on m ass, which is an important discriminant of 



galaxy formation scenarios (jPeletier et al, 



19901). 



In the present paper, we derive the PHP in optical and NIR wavebands for a large sample 
of ETGs belonging to the cluster of galaxies GIG 1008-1224, also known as M S1008-1224 



(hereafter MS 1008), a rich X-ray bright cluster at z = 0.306 (iLewis et al. 
detected in the Einstein Medium Sensitivity Survey ( Gioia fc LuppinJll994 ). For this cluster. 



199^, originally 



a unique wealth of multi- wavelength data is available from the ESQ archive ^, including 
very deep UBVRIJHK photometry taken with VLT FORS and ISAAG and NTT SOFI 
instruments in excellent seeing conditions. The data used for the present analysis consist of 
the R-, I- and K-band photometry of MS 1008. This unique data-set allows us to perform 
for the first time a homogeneous and accurate multi-wavelength analysis of the PHP at 
intermediate redshifts. In the present paper we study (i) the characteristics of the PHP 
relation (i.e. its coefficients, scatter and shape) at z ~ 0.3, accounting for selection effects 
as well as for the fitting procedure; (ii) the waveband dependence of the PHP; and (iii) 
the variation of the PHP coefficients from z ~ 0.3 to z ~ 0. We consider the following 
representation of the PHP: 



http: / / archive.eso.org 
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logTe = alogn + 6 </i>e +c, (1) 

where a and b are the slopes and c is the zero-point of the relation. 

The layout of the paper is as follows. In Section 2 we describe the samples used in 
the analysis, while in Section 3 the surface photometry is presented. Section 4 deals with 
the problems related to the fit of the PHP. In Section 5 we derive the PHP at z ~ 0.3 
in both optical and NIR wavebands. In Section 6 we present the waveband and redshift 
dependence of the PHP. The variation of the PHP coefficients with redshift is analyzed 
trough a comparative analysis of the PHP at 2; ~ 0.3 with that of ETGs in clusters at z ~ 0. 
The main results are then discussed in Section 7. A summary follows in Section 8. In the 
present paper we assume a ACDM cosmology with Hq = 70 Km s^^ Mpc^\ Q^n = 0.3 and 

= 0.7). 



2 THE SAMPLES 



The galaxies used for the present study belong to the cluster MS 1008 ai z = 0.306 and 
were selected on the basis of a large wavelength baseline, including UBVRIJHK photome- 



try for a field o 



see 



Gioia et al. 



6.8' X 6.8' around the cluster center (aj2ooo=10:10:34.1, 5j2ooo=-12:39:48, 



19901 ). This dataset, which was retrieved from the ESO archive and reduced 



by the authors as described in Covone et al. (2005, in preparation, hereafter CLB05), was 
used to estimate photometric redshifts for all the sources in the cluster field. Details on pho- 
tometric redshifts can be found in CLB05, while we give here only the relevant information. 
In CLB05 we compared the photometric redshifts with spectro scopic redshifts available for 
~ 70 galaxies in the cluster field from the CNOC survey ()Yee et al.Ml998h . The mean 



offset and the rms of differences between spectroscopic and photometric redshifts turned 
out to be 0.004 ± 0.005 and (fe~0.04 respectively. The quantity 5z gives an estimate of the 
typical accuracy on photometric redshifts. We selected as cluster members the objects with 
photometric redshifts in a range of ±0.08 (i.e., ±2(fe) around the cluster redshift. 

For the present study, we use only the i?-, /- and A"— band photometry of MS 1008, 
whose depth and resolution are suitable to obtain structural parameters for a fair sample 
of galaxies. The R- and J-band data were taken with VLT-F0RS2 at ESO, while the K- 
band imaging, including ESO-VLT and ESO-NTT photometry, was taken with ISAAC 
and SOFI instruments respectively. The R- and /- band images cover a field of 6.8' x 6.8' 
with a pixel scale of 0.201 "/pixel, while the A'-band SOFI, hereafter K{S), and A'-band 
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ISAAC, hereafter K{I), images cover two fields of 5' x 5' and 2.5' x 2.5' around the cluster 
center with pixel scales of 0.288 "/pixel and 0.147 "/pixel respectively. The exposure times 
of the R-, /-, K{S)- and K{I)- images were 1.5, 1.125, 5.7 and 1.4 hours, while mean 
seeing sizes were 0.8", 0.9", 0.8" and 0.5" respectively. We derive the PHP for four different 
samples of galaxies, in the R-, /-, K{S)- and i^'(J)-band images respectively. Galaxies in 
each sample were selected as follows. In order to obtain reliable structural parameters, 
we selected galaxies down to a given magnitude cut, which was established by deriving 



La Barbera et al 



2mm 



structural parameters from simulated galaxy images (see sec. 3. 3 of 
for details). For each waveband, we estimated the magnitude limit to which systematic 
uncertainties on effective radius and Sersic index are expected to be negligible. Hence, we 
considered only cluster members brighter than R = 22.0, / = 21.1, and K = 17.8 in the 
R-, I-, and fT-band images respectively. We also excluded those objects very close to bright 
saturated stars in the field, since structural parameters of those galaxies would have large 
systematic uncertainties due to background subtraction problems. A further selection was 
applied a posteriori, by selecting as ETGs the objects with Sersic index n>2 and excluding 
objects with uncertainties on structural parameters greater than 100%. The adopted Sersic 
inde x cut corresponds to ex clude disk-dominated systems from the present analysis (see 



e.g. 



van Dokkum et al. 



1998^ . The fraction of excluded objects was smaller than 7% for 



each band. The final R, I, K{S) and K{I) samples include = 129, A^ = 123, A^ = 68 
and A^ = 50 galaxies respectively, with A^ = 112 objects in common to the R— and /—band 
images, and A^ = 38 galaxies are in common to the K{S) and A'(J)samples. The position 
of these galaxies in the cluster field is shown in Fig. ^ 

We point out that the above selections allow us to obtain for each sample two sharp 
selection cuts in the space of structural parameters, i.e. the magnitude limit and the Sersic 
index cut. This is crucial for an accurate correction of PHP coefficients for selection effects 
(Section 112 



3 DERIVATION OF STRUCTURAL PARAMETERS 



Structural parameters were derived for galaxies in the R-, /-, K (S)- and A"f/)-band ini ages 



La Barbera et al 



2002 and 



of MS 1008 by using the two-dimensional fitting method (see 
references therein). The surface brightness of galaxies was modeled by the Sersic law (|Sersic 
19681 ): 
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Figure 1. FORSl R-band image of MS 1008. The image covers an area of ~ 6.8' X 6.8' around the cluster center. Galaxies 
included in the optical (R and I bands) and in the K-band samples are marked by circles and squares respectively. 



yu(r) = fiQ + 2.5bn log e ■ (r/ r^) 



l/n 



(2) 



where bn is a constant, defined in such a way that r^. is the effective (half-hght) radius of the 
galaxy, r is the equivalent radius, n is the Sersic index (shape parameter), and is the central 



surfa ce brightness. The constant 6„ can be estim ated bv a power- law in n fICiotti fc Bertin 



19991 1 with 6„ ~ 2n — 1/3 in first approximation ()Capacciolilll989l ) . or, with higher accuracy 
(< 1%), by the function 6„ ~ exp[(0.6950 + lnn) -0.1789/n] (see LG M99). The mean surfa ce 
brightness within rg, </i>e, is given by the following formula fe.g. lCiotti fc Bertinlll999| ): 



8 F. La Barbara et al. 



= -2.51og(27r) +;Uo - 51og(re) + 

+5n log 6„ - 2.5 log r(2n) 
= -2.5 log(27r) - 5 log(re)+ < /i >e, 



x' = E[Gi,j-(bg + S®PSFi,)f , (4) 



where mx is the total magnitude of the galaxy and F is the complete gamma function. 
Galaxy images were fitted with seeing-convolved Sersic models, minimizing the function: 

,2 _ ^ rr^ n i a -dot? m2 

where the symbol ® denotes convolution, Gij is the galaxy surface brightness at pixel 
bg is the local background value, while S and PSF are the Sersic and the PSF models 
respectively. Seven output parameters were provided from the fitting process: the center 
coordinates, the effective radius, r^,, the mean surface brightness, </i>e, the Sersic index, n, 
the position angle, PA, of major axis and the axis ratio, b/a. For each galaxy, neighbor ob- 
jects were masked automatically using the ellipticity, pos ition angle and isoarea parameters 



provided by running SExtractor (jBertin fc Arnoutlll996^ on the corresponding images. In 



the cases where SExtractor failed to provide reliable estimates of these parameters, such as 
in crowded regions and/or in the neighborhood of extended sources, masking was performed 
interactively. Very close galaxies were fitted simultaneously. 

The PSFs of the R, I, K{S) and K{I) images were modeled by a sum of 2D Moffat 
functions, taking into account deviations of stellar isophotes from circular symmetry. Details 
on PSF modeling are given in Appendix 1X1 Uncertainties on log re, <yU>e and logn were 
estimated by using numerical simulations and by comparing structural parameters among the 
different wavebands. Details can be found in AppendixiBl Structural parameters^ for galaxies 
in the R-, I- and i^'-band samples are reported in Table 1. Since structural parameters turned 
out to be fully consistent between the K{S) and K{I) samples (see Appendix EJ, in Table 1 
we report the averages of the K{S) and K{I) structural parameters, which were computed 
by weighting each value with the inverse square of the corresponding uncertainty. 



4 FITTING THE PHOTOMETRIC PLANE 

We derived the coefficients of the PHP by accounting for (i) selection effects, i.e. the 
cuts in Sersic index and magnitude; and (ii) the correlated uncertainties on structural 

^ Table 1 is only available in electronic form or via anonymous ftp to cdsarc.u-strasbg.fr at the CDS. 
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Table 2. Acronyms used for fitting methods. Since the fit coefficients were corrected for selection effects (see Section 14.21 . a 
prefix C was used for each acronym. 



CLSiogT-e Corrected weighted least square fit with dependent variable logrc. 

CLS<^^ Corrected weighted least square fit with dependent variable <fi>c 

CLSiog n Corrected weighted least square fit with dependent variable log n 

CORLS Corrected orthogonal weighted least-square fit 

CBLS Corrected bisector least square fit 

CAMLS Corrected arithmetic mean of the least-square coefficients 

CGMLS Corrected geometric mean of the least-square coefficients 



parameters. Since different fitting meth ods are not equally affected by these effects (see 



La Barbera. Busarello. Capacciolil 1200(1 LBCOO), the PHP coefficients were obtained for 
each sample by using different fitting methods, as described in Section 14.11 Correction for 
selection effects was performed by a Monte Carlo technique as detailed in Section 14.21 



4.1 Regression methods 

The PHP coefficients were derived by using seven different fitting procedures: three weighted 
ordinary least-square fits, adopting as dependent variable one of the quantities log re, <fi>e 
and logn; the orthogonal weighted least square fit, where the root mean square (rms) of 
residuals perpendicular to the plane are minimized; the bisector least square (BLS) fit (see 
LBCOO); the arithmetic and the geometric means of the ordinary least-square coefficients. In 
the following, acronyms will be used for different fitting methods as summarized in Table El 
In each fitting procedure, we derive the coefficients a, b and c of the plane (see Eq. and 
its intrinsic dispersion a\ For the ordinary and orthogonal least square fits, these quantities 
were derived by minimizing the following expression: 

= -lnL = 



i:^;;(3SwF + V21nh(a,6,c)^ + (aO^; 



(5) 



where L is the likelihood function, are the residuals about the plane, are the uncer- 
tainties on Tfc, and is the intrinsic dispersion of the PHP along the direction, dmin? where 
residuals are minimized. The terms are the uncertainties on and were obtained by pro- 
jecting the covariance matrix of measurement errors on structural parameters along dmin- 
We point out that although this procedure allows each point to be weighted with the corre- 
sponding uncertainties, it does not correct exactly for biases on fitting procedures which are 
due to the correlation among such uncertainties. It is always possible, in fact, to vary the 
covariance matrix of measurement errors on log re, <yU>e and logn, without changing its 
projection on a given direction dmin. In order to estimate the amount of bias which is due 
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to the correlation of the uncertainties, we adopted the MIST fits (LBCOO) which provides 
unbiased values of ordinary least-square coefficients by using the mean covariance matrix of 
uncertainties on the three variables. The bias turned out to be always smaller than 5% for 
each sample, and therefore we did not apply corrections for this effect. In order to derive the 
quantities a, 6, c and a', we minimized Eq. Elby using a Levenberg-Marquardt algorithm. For 
what concerns the CBLS, CAMLS and CGMLS fits, the PHP coefficients were calculated 
from those of the CLS fits. For this reason, a different weight is given to each galaxy also 
in these fitting methods. Uncertainties on a, 6, c and o"^ were estimated by the bootstrap 
method, applying = 2500 bootstrap iterations. 



4.2 Selection effects 

Selection effects were corrected for by using numerical simulations, producing distributions 
of points in the space of structural parameters which resembled those of real galaxies. Points 
in each simulation were generated as follows. Magnitudes were assigned according to the lu- 
minosity function, which was modeled as a Schechter function with slope a and chara cteristic 



magnitudes M* . The val ues of a and M* were drawn from 



R and / bands, and from 



de Propris et a1 



Busarello et al. 



(2002) for the 



( 19991 ) for the K band. Effective radii and mean 



surface brightnesses were obtained according to the luminosity-size relation, whose slopes, 
zero-point and scatter were derived from the data of MS 1008. Sersic indices were assigned 
using the PHP coefficients, a, h and c, and its logn scatter, criog„. The values of a, 6, c and 
c"iogn were chosen by an iterative procedure. For each iteration, a simulated distribution in 
the space of structural parameters was constructed by imposing the same selection effects, 
i.e. the magnitude and Sersic index cuts, of the real samples, and the PHP coefficients were 
computed by using all the fitting methods described in Section 14.11 The values of a, 6, c 
and (TiognWere modified until the coefficients derived from the simulated samples matched 
those of the real samples for all the fitting methods. This procedure allowed us to achieve 
an excellent match, with an accuracy better than ~ 5% for all the fitting coefficients. Once 
the values of a, 6, c and aiogn were chosen, we estimated the relative variations of PHP co- 
efficients of all the regression methods after selection cuts were removed from the simulated 
samples. The variations were used as correction factors for the PHP coefficients of MS 1008. 
As an example, in TableElwe report the correction factors for the i?-band sample of MS 1008. 
We see that the bias can be as large as ~ —35% for the coefficient a of the CLSwn fit, and 
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Table 3. Bias on the /?-band PHP cocfBcients due to selection effects. 





Sa 


Sb 


5c 


(5a-iog rc 






a 


b 


c 


CTlog re 


log rc 


CLSlog ro 


0.22 


-0.07 


-0.02 


0.12 


0.08 




0.33 


0.08 


0.09 


0.18 


0.16 


CLSlog n 


-0.35 


0.08 


-0.01 


-0.16 


-0.03 


CORLS 


-0.10 


-0.01 


-0.03 


0.02 


0.09 


CBLS 


0.10 


0.00 


0.01 


0.11 


0.12 


CAMLS 


0.02 


0.03 


0.03 


0.09 


0.12 


CGMLS 


0.19 


0.03 


0.05 


0.11 


0.08 



Table 4. Coefficients of the PHP in the _R band. 



log i\ 



CLSiogre 1.00 ±0.09 0.184 ±0.010 -4.70 ±0.21 0.17 ±0.01 0.15 ±0.01 

CLS^>„ 0.90 ±0.14 0.303 ±0.019 -7.07 ±0.38 0.21 ± 0.02 0.19 ±0.03 

CLSiogn 1.40 ±0.16 0.157 ±0.025 -4.22 ± 0.40 0.18 ±0.02 0.16 ±0.02 

CORLS 1.10 ±0.11 0.175 ±0.012 -4.48 ± 0.22 0.16 ±0.01 0.14 ±0.01 

CBLS 1.07 ±0.06 0.219 ±0.009 -5.41 ±0.19 0.17 ±0.01 0.14 ±0.01 

CAMLS 1.09 ±0.08 0.215 ± 0.010 -5.34 ±0.20 0.17 ±0.01 0.15 ±0.01 

CGMLS 1.07 ±0.09 0.201 ± 0.017 -5.04 ±0.34 0.17 ±0.01 0.15 ±0.01 



that it strongly depends on the fitting method, being smaller than ~ 10% in absolute value 
for the CORLS, CBLS and CAMLS fits. The uncertainties on the correction factors were 
estimated by changing the input parameters of the simulation algorithm (e.g. the value of 
a) according to their uncertainties, and repeating the iteration procedure. We found that 
the bias on the PHP coefficients varies at most by ~ 4% for the coefficient a in the CLS fits. 
This is smaller with respect to the typical uncertainties on PHP coefficients, and was there- 
fore neglected. We point out that the correction procedure was repeated independently for 
each waveband of MS 1008, without making any assumption a priori on how the properties 
of galaxies in the space of structural parameters can vary among the different wavebands. 



5 PHOTOMETRIC PLANES OF MS 1008 

The distribution of galaxies in the space of structural parameters is shown in Fig. |21 where we 
show the distribution of galaxies in the logre-</i>e, logn-</i>e and logn-logre planes, 
and a 3D view of the PHP, with the corresponding magnitude and Sersic index cuts repre- 
sented by two orthogonal planes. An edge projection of the plane is also shown in Fig. 01 
We see that galaxies follow a well-defined PHP at z = 0.3, with Sersic indices that increase 
towards lower surface brightness values and larger effective radii. 

The coefficients of the PHP in the R, I and K bands for the different fitting procedures 
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Figure 2. PHP of MS 1008 in the R-band. The upper-left panel shows a 3D view of the plane, where black circles mark 
positions of galaxies (N = 129) in the space of structural parameters. The shaded surface is the PHP, as defined by the CBLS 
fitting method, while the vertical and horizontal blue grids depict the selection cuts in magnitude and Sersic index respectively. 
The shading is realized in such a way that color intensities are inversely proportional to logn. Note that the plot of the PHP 
is transparent, and therefore black points below the surface appear as partly obscured symbols. The upper— right and the lower 
panels show the logrc-</t>c, logri-</^>c and logn-logre projections of the PHP. 



are reported in Tables|U[Sland|nirespectively, together with the estimate ciJog rc intrinsic 
dispersion in logrg about the plane. All these values were corrected for selection effects as 
detailed in Section 14.21 The i^-band coefficients were obtained by combining the values 
obtained for the i^(S')and K{I) bands, as discussed below (Section I5.2|l . For the CLS and 
CORLS fits, the quantity crf^g,,.^ was estimated by projecting in Eq. El along the direction 
logrc, while for the other fitting procedures we subtracted in quadrature to the rms of logrg 
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residuals the amount of scatter due to measurement errors. This was performed by taking 
into account the correlation of uncertaintities on logrg, </i>e and logn. 

By looking at Tables EJ El and El we see that (i) the coefficients of the PHP obtained by the 
various fitting procedures are significantly different, and (ii), whatever the regression method 
is, the PHP has significant intrinsic dispersion. As discussed by LBCOO for the FP relation, 
the existence of intrinsic scatter for a bivariate relation and our ignorance on its origin imply 
that different fitting methods do not necessarily provide consistent results. We note, however, 
that the dependence of PHP coefficients on the fitting procedure is particularly significant 
only for the CLS regressions, which give a special role to the variable whose residuals are 
minimized during the fit. On the other hand, treating equally all the variables, as in the 
CORLS, CBLS, CAMLS and CGMLS fits, gives much more robust estimates, allowing us to 
obtain a stable determination of PHP coefficients. We also note that the CLS<^^ regression 
provides a lower value of a and a higher value of b with respect to the CLSiogrc method. 
On the contrary, the CLSiogn regression produces a higher value of a and a lower value of b 
with respect to the CLSiogrc fit- Since the CBLS, CAMLS and CGMLS fits are based on an 
average of the CLS coefficients, the corresponding coefficients are very close to those of the 
CLSiogrc method. Since the coefficients of the bisector fit have a smaller relative uncertainties 
with respect to the other fitting methods, we will refer to these coefficients in the following 
analysis f Section 16. 2|) . It is also interesting to note that the correction for selection effects 
is negative for the CORLS fit while it is positive for the CBLS, CAMLS and CGMLS 
regressions, and that, therefore, if selection effects would have not taken into account, the 
difference among the CORLS and the CBLS, CAMLS and CGMLS fits would have been 
significant. 

5.1 The PHP in optical wavebands 

A proper comparison of PHP coefficients between different samples has to account for the fact 
that their uncertainties are correlated. For this reason, we chose to compare simultaneously 
the slopes a and b of the PHP and one slope, a, with its zeropoint c. The comparison of 
b versus c as well as other plots with pairs of quantities from Tables IH and El do not add 
further information to the discussion and are not shown in the following. In Fig. |31 we plot 
the R- and /- band coefficients. Taking into account the mean galaxy colour, i?— / ~ 0.75 at 
z ~ 0.3, we expect an offset between the I- and R- band PHP zero-points of A^^j ~ b-R—I. 
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a-logn + b-</i,>g + c 

Figure 3. Edge-on view of the H-band PHP of MS 1008. The plot shows the A'' = 129 galaxies of the i?,-band sample. Ellipses 
denote uncertainties on structural parameters and correspond to la confidence contours. 

In order to remove this effect from tlie comparison, we subtracted tlie quantity A^_/ from 
tlie I- band values of c wliicli are sliown in Fig. HJ We note tliat since tlie value of A^_/ is 
proportional to 6 ~ 0.2, uncertainties on R—I do not affect the comparison in Fig. |3J 

As shown in the figure, the PHP coefficients are fully consistent between the R and / 
bands. This is clearly in agreement with the fact that structural parameters do not show 
significant variations among optical wavebands (see Appendix^, and is due to the fact that 
the R and / bands at z = 0.3 sample a very similar spectral region, where differences in 
stellar population properties of galaxies, such as their internal colour gradients and the colour 
magnitude relation, are negligible. For what concerns the dispersion around the plane, we 
see from Tables 0] and El that this is also consistent between the R and I bands, amounting on 
average to ~ 0.18 dex (~ 41%) in logrg. The same result holds for the intrinsic dispersion of 
the PHP, which amounts to ~ 0.14 — 0.15 dex (~ 32 — 35%) in logrg. We note that although 
for the CLS<^^ and CLSiog„ fits the value of crlogr^ is slightly larger with respect to that of 
the other fitting methods, the corresponding uncertainties of cJogre larger, making this 
difference not particularly significant. We also note that only a few percent of the observed 
dispersion about the PHP are explained by the the measurement errors on logr^, <fi>e 
and logra. As shown in Fig. IHl in fact, uncertainties on structural parameters are strongly 
correlated in a direction which is almost parallel to the PHP. 

We also analyzed the presence of a possible curvature in the distribution of ETGs in the 
space of logTe, </i>e and logn, by considering correlations among log re residuals about 
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c 

Figure 4. Comparison of tlie optical PHP coefficients. The R- and /- band coefficients are marked by asterisks and crosses 
respectively, while the fitting procedure are shown in different colors, as indicated in the upper panel. The solid (dashed) ellipses 
mark la confidence contours for each pair of R(I)-band coefficients (see text). The dashed lines in the upper and lower panels 
show the direction of the mean correlation of uncertainties on a and b, and b and c respectively. 



the plane and each of the three variables. We found that such correlations are very weak, 
the Spearman's rank coefficients amounting to 0.14 for the correlation between residuals 
and logTe, and to -0.2 for the other correlations with <fi>c and logra. These coefficients 
can be easily explained by the selection cuts in the space of structural parameters. We also 
fitted the -R-band PHP by using log as dependent variable and by adding to Eq. ^ one 
among all the possible quadratic terms which can be constructed from the two quantities 
</i>c and log 72. All these terms were found to be consistent with zero at the level of lAa, 
the log Tg scatter of the PHP decreasing by less than 1% in each case. Our data, therefore, 
indicate that there is no significant departure of the PHP from a flat relation. 
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Figure 5. The same of Fig. HI for tlie K{S) and K{I) coefficients. Asterisks and crosses mark the K{S) and K{I) PHP 
coefficients respectively. EUipses denote l.Scr confidence contours. 



5.2 The Ji-band plane 

In order to derive the PHP of MS 1008 in the K band we apphed separately the different 
fitting procedures to the K{S) and K{I) structural parameters. The comparison of the 
K{S) and K{I) coefficients is shown in Fig. |3[ where the same quantities as Fig. |^ are 
shown. The important outcome is that the PHP coefficients of the NIR samples are fully 
consistent, with a confidence level < l.Scr. The same result holds for the dispersion around 
the plane. This is a reassuring result since the K{S) and K{I) parameters were obtained 
from images with different resolution and seeing. 

Since the coefficients and dispersions of the K{S) and K{I) PHPs turned out to be 
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Table 5. Coefficients of the PHP in the I band. 



flog CTjQg 



CLSiogr„ 0.90 ±0.15 0.194 ±0.013 -4.73 ± 0.25 0.18 ±0.01 0.14 ±0.01 

CLS^>^ 0.75 ±0.20 0.331 ±0.026 -7.32 ± 0.51 0.22 ± 0.02 0.20 ±0.02 

CLSiogn' 1.66 ±0.24 0.124 ±0.029 -3.70 ±0.46 0.25 ± 0.04 0.10 ±0.05 

CORLS 1.19 ±0.13 0.182 ±0.015 -4.53 ± 0.26 0.19 ±0.02 0.14 ±0.01 

CBLS 1.03 ±0.09 0.232 ±0.013 -5.52 ± 0.24 0.19 ±0.01 0.15 ±0.01 

CAMLS 1.17 ±0.12 0.218 ±0.014 -5.33 ± 0.24 0.20 ±0.01 0.16 ±0.01 

CGMLS 0.95 ±0.16 0.186 ±0.025 -4.55 ± 0.51 0.19 ±0.02 0.15 ±0.02 



Table 6. Coefficients of the PHP in the K band. 



flog i-o flog 7\ 



CLSiogre 0.83 ±0.13 0.209 ±0.010 -4.51 ±0.19 0.18 ±0.01 0.15 ±0.02 

CLS^>„ 0.75 ±0.17 0.351 ±0.020 -6.83 ± 0.32 0.21 ± 0.02 0.17 ±0.02 

CLSiogn 1.39 ±0.29 0.151 ±0.033 -3.51 ± 0.39 0.27 ±0.05 0.23 ± 0.05 

CORLS 0.90 ±0.14 0.192 ±0.010 -4.18 ±0.19 0.19 ±0.01 0.17 ±0.02 

CBLS 0.97 ±0.08 0.245 ± 0.010 -5.13 ±0.17 0.20 ±0.01 0.15 ±0.02 

CAMLS 1.00 ±0.15 0.235 ±0.013 -4.98 ±0.19 0.21 ± 0.02 0.19 ±0.02 

CGMLS 0.91 ±0.13 0.212 ± 0.027 -4.57 ±0.49 0.19 ±0.02 0.16 ±0.02 



consistent, we combined their values by a weighted mean. The resulting coefficients are 
shown in Table IHl 



6 PHP DEPENDENCE ON WAVEBAND AND REDSHIFT 
6.1 The PHP in optical and NIR wavebands 

The R- and K- band coefficients of the PHP are shown in Fig. IHl In order to remove the 
mean galaxy colour from the comparison, we proceeded as in Section 15.11 by subtracting 
from the coefficient c in the K band the term Aji^k = —b ■ R—K, with R — K 3.2 at 
z ~ 0.3. Since the /-band coefficients are consistent with those in the R band, they do not 
add further information to the discussion and are not shown in the figure. By looking at 
Fig. inland at TablesEJElandini we see that, whatever regression method is adopted, the PHP 
coefficients turn out to he fully consistent among the optical and NIR wavebands. We note 
that since the Sersic index and the magnitude cuts affect each fitting procedure in a different 
way, the independence of the above result on the ffiting method makes it very robust with 
respect to selection effects in the samples. Considering the CORLS, CBLS, CAMLS and 
CGMLS ffis, we find a ~ 1 and h ~ 0.2 both in the optical and NIR wavebands. We also 
note that the intrinsic scatter about the K-band plane is fully consistent with the value of 
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Figure 6. Comparison of the PHP coefficients for the R and K bands (asterisks and crosses respectively). SoHd (dashed) 
eUipses denote Icr confidence contours for R(K)-band coefficients. 

~ 0.14 — 0.15 dex found for the optical PHP, provided that measurement uncertainties are 
taken into account. 

6.2 The PHP at 2 ~ 0.3 and z ~ 

In order to address the redshift dependence of the PHP, we compared our resuhs with 
those of GRA02, who derived the PHP at 2; ~ by using 5-band structural parameters for 
= 38 ETGs in the Virgo and Fornax clusters. We note that the R band aX z = 0.3 is 
closed to i?-band restframe, and therefore the data of GRA02 and those of MS 1008 cover 
approximately the same restframe wavelengths. To perform a homogeneous comparison, we 
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Figure 7. Comparison of the i?-band PHP slopes of MS 1008 (crosses) with those of nearby galaxies in the B band from 
GRA02 (asterisks). Ellipses denote l.Scr confidence contours. 

used the data in table 1 of GRA02 and re-derived the PHP coefficients using the regression 
methods described in Section I^^Tl Selection effects were taken into account by constructing 
PHP simulations as described in Section with the same magnitude and Sersic index cuts 
as the GRA02 sample, i.e. B = 15 and n = 1.5, respectively. We note that, since the GRA02 
sample is not complete in any sense (see sec. 2 of GRA02), the above procedure provides 
only a rough correction of the selection biases. Because of the wide magnitude range of 
the GRA02 sample (~ 5.3 mag), the corrections turned out to be quite small, amounting 
to 5% on average, and to at most 10% for the coefficient a in the CLSiogn fit. The bias 
corrected PHP coefficients at 2; ~ are given in Table [7| The intrinsic dispersion of the 
PHP was computed assuming a typical uncertainty of ~ 0.1 dex in logn (see sec. 2 of 
GRA02), and adopting as mean covariance matrix of uncertainties on structural parameters 
that of galaxies in MS 1008, which was re-scaled to match the uncertainty on logn at ^ ~ 0. 
The slopes of the PHP obtained by the CORLS, CBLS, CAMLS and CGMLS ffis can be 
compared with the values of a = 0.89 ± 0.14 and b = 0.24 ± 0.036, that were derived from 
GRA02 by the bisector ffi, treating equally all the three variables. The values in Table [7| 
are fully consistent with those of GRA02. The PHP slopes of the GRA02 sample and those 
of MS 1008 are compared in Fig. [3 where we see that the values of a and b at 2; ~ are 
fully consistent with those at z ~ 0.3. We note that due to the smaller sample size at 2; ~ 
the corresponding uncertainties on a and b are larger, particularly for some of the ffiting 
procedures (e.g. the CLS<^p fit). As shown in Table [71 the intrinsic dispersion around the 
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Table 7. B-band coefficients of the PHP for the GRA02 sample. 



CLSiogr„ 0.97 ±0.16 0.221 ±0.039 -4.87 ±0.80 0.17 ±0.01 0.14 ±0.02 

CLS^>^ 0.58 ±0.27 0.433 ±0.078 -9.00 ± 1.55 0.20 ±0.04 0.17 ±0.04 

CLSiogn 1.31 ±0.24 0.135 ±0.052 -3.18 ±0.98 0.21 ± 0.04 0.15 ±0.03 

CORLS 1.14 ±0.15 0.180 ±0.031 -4.07 ±0.64 0.19 ±0.02 0.14 ±0.01 

CBLS 0.92 ±0.12 0.279 ± 0.036 -5.92 ± 0.71 0.17 ±0.01 0.13 ±0.01 

CAMLS 0.97 ±0.12 0.263 ± 0.034 -5.58 ± 0.66 0.17 ±0.01 0.13 ±0.01 

CGMLS 0.85 ±0.30 0.212 ± 0.060 -4.45 ± 1.25 0.18 ±0.05 0.16 ±0.06 

plane does not change significantly, provided that measurement uncertainties are taken into 
account. 

Taking advantage of the fact that the PHP slopes are consistent between 2; ~ and 
z ~ 0.3, we re-computed the coefficient c at 2; ~ and z ~ 0.3 by fixing the slopes of 
the PHP to the values of a and h obtained for MS 1008 by the CBLS fit, i.e. a = 1.07 
and b = 0.219. The CBLS fit was chosen because of the smaller uncertainties of bisector 
coefficients with respect to those of other fitting procedures. The value of c at 2; ~ was 
also corrected for magnitude and Sersic index cuts as described in Section 14.21 Due to 
the robustness of the CBLS regression with respect to selection effects (see Table EI), the 
correction turned out to be very small (~ 0.2%). The coefficient c amounts to — 4.737±0.028 
at ;z ~ and -5.337 ± 0.014 for MS 1008, differing by A(c) = -0.60 ± 0.03. The value of 
A(c) can be analyzed by taking the average of the difference between the PHP equations 
(Eq. HI) at z ~ 0.3 and z ~ 0: 



A(c) = A(logre) - a ■ A(logn) - b ■ A(</i>e). (6) 

Assuming that, on average, structural parameters vary from 2 ~ to z ~ 0.3 only because 
of the luminosity evolution of stellar populations, we can re-write the previous equation by 



setting A(logre) = and A(logn) = (see discussion in Section mj) : 

1 + 



A(c) = -log(dA)-6 
where the term lOA 



lOA 



log 



0.3) 



(7) 



1.008 

log (y^) accounts for the surface brightness dimming between the 
redshift of MS 1008 and that of the GRA02 sample^, dp^ is the angular diameter distance^ 
corresponding to z = 0.306, while (-B0.0 ~ R0.3) is the difference between the B magni- 
tude at 2; ~ and the -R-band magnitude at 2; ~ 0.3. Eq. [7| was used to compare the 

^ For the GRA02 sample of Virgo and Fornax galaxies, we assumed z = 0.008. 

^ The term — log(dA) accounts for the fact that effective radii of galaxies in MS 1008 are given in arcsec, while for the GRA02 
sample Tc is given in kpc. 
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Figure 8. The difference of the PHP zero-point, A(c), between z = 0.3 and z = is compared with predictions of stellar 
population models. The dash-dotted line and the gray hatched region indicate the mean value of A(c) and the corresponding 
±1.5cr confidence interval. The dotted lines mark the intersections of the two models with the hatched area. 

meas ured value of Ajc) with p redictions of stellar population models from the G1SSEL03 



code ( Bruzual fc Chariot 



20031 ). We considere d two different stellar population models with 
J II 1 

solar metallicity and a Scalo IMF ()Scaldll986l ). the first one being a simple stellar population 
(SSP), the second one having an exponentially declining star formation rate with time scale 
r = 1 Gyr. For both models, we computed the value of A(c) from Eq. [7|by deriving the 
B- and R- band magnitudes at ^ = and z = 0.306 respectively, for different formation 
redshifts Z{. In Fig. |Hl we show the values of A(c) for the two stellar populations as a func- 
tion of Zf. The measured value of A(c) is marked by the dash-dotted line in the plot, while 
the hatched area is the confidence interval defined by measurement uncertainties on A(c). 
The figure shows that both models cross the hatched region and are consistent, therefore, 
with the measured value of A(c). The crossing points define lower limits for the formation 
redshift of galaxy stellar populations. For the SSP and r = 1 Gyr models, we obtain Zf > 1 
and Zf > 2 respectively. 



7 DISCUSSION 

We have shown that cluster ETGs at 2; ~ 0.3 follow a tight correlation among log re, </i>e 
and logn, with an intrinsic dispersion of ~ 32% in r^. Our data indicate that the PHP has 
no significant curvature. In Section 17.11 we discuss whether stellar populations can be the 
origin of the intrinsic dispersion about the plane, while in Section 17.21 we attempt to find a 
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possible explanation for the absence of curvature in the PHP. Sections 17.31 and 17.41 deal with 
the waveband and redshift dependence of the PHP coefficients. 



7.1 The intrinsic dispersion of the PHP 

In Section El we showed that the PHP relation has significant intrinsic dispersion, amounting 
to (T}og^^ ~ 0.14 dex (~ 32%) both in optical and NIR wavebands. This dispersion is fully 
consistent with that found by GRA02 for ETGs in nearby clusters. The intrinsic scatter 
of the PHP turns out, therefore, to be larger with respect to that of the FP (see values 
reported in Section although it is s maller with respect to that of other co rrelations 
among galaxy parameters. For example, iPahre. Dioreovski fc de Carvalhd (|l998l ) analyzed 
different correlations among photometric and spectroscopic properties of nearby ETGs. They 
found that the Kormendy relation (KR) and the Mg2 FP have a log re scatter of ~ 0.23 dex 
and ~ 0.17 dex respectively, where the Mg2 FP relation was constructed by replacing the 
velocity dispersion term in the FP with the Mg2 line-strength. It is interesting to compare 
these values with those we can obtain for MS 1008. A KR fit to the -R-band sample of MS 1008 
gives a log re dispersion of ~ 0.2 dex, that is ~ 14% higher than the scatter of the PHP. This 
result is consistent with what found by LMB04. Since we do not have Mg2 line-strengths, 
we constructed a 'photometric' FP bv replacing y elocitv dispersions with optical-NIR galaxy 
colours (see e.g. Ide Carvalho fc Djorgovskilll989l ). Using I—K colours available for = 105 
galaxies from the i?-band sample of MS 1008, we obtain logrg oc (0.25 ± 0.02)- < /i >e 
+ (0.23 ±0.07) • (I—K), with a log re rms of 0.175 dex, which is ~ 8% larger with respect to 
the PHP, but significantly smaller with respect to the KR. Since effective parameters turn 
out to correlate with both colours and Sersic indices, we also tried to construct a photometric 
hyperplane with log re, </i>e, logn and I~K. We found, however, that the I — K term 
in such a hyperplane is only marginally significant (1.5cr), and that the scatter about the 
relation decreases only by 1% with respect to that of the PHP. This result is also shown in 
Fig. El where we plot log re residuals about the PHP versus residuals from the I — K vs. K 
CM relation of MS 1008. Such 6 — 6 diagrams have be en already applied in sev eral works 
to investigate the origin of residuals about the FP (e.g. iPrugniel fc SimienI 119961 ). The plot 
clearly shows that there is no correlation among residuals to the PHP and CM relations, 
the corresponding Spearman's rank coefficient amounting to —0.12 ±0.09. This implies that 
either the scatter of the PHP does not origin from stellar population parameters or that a 
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Figure 9. Relation among residuals to the CM and PHP relations. 

complex combined effect of such parameters (e.g. age and metallicity) is acting in such a way 
that no correlation appears in the 6-S diagram. This issue could be addressed by correlating 
PHP residuals with line-strength indices. 



7.2 Why a plane with logrg, <yU>e and logra? 

As mentioned in Section^ MLCOO and MLCOl found that ETGs populate a curved manifold 
in the space of structural parameters, showing that the existence of such a PHP can be 
explained by the specific entropy of ETGs, s, being an increasing function of their mass, M. 
MLCOl also showed that the existence of a plane in the space of logrg, /xq and logra can 
be explained by the limited range considered for logn in previous works (with logn > 0.15, 
see fig. 10 of MLCOl). Since the variables used by MLCOl as well as the quantities logrg, /io 
and logn are not linearly related to logrg, </i>o and logn (see Section [T]), an intriguing 
question arises on what the origin is of the flatness of the PHP in the logrg, </i>g and 
logn space. 

To address this issue, we tried to understand if there is some simple relation between 
the specific entropy of ETGs and the three quantities logrg, </i>g and logn. This was 
achieved by considering the same galaxy models described in MLCOO and MLCOl, that is 
spherical, non-rotating systems, with negligible radial gradients of the mass-to-light ratio 
(i.e. one component models). We point out that attempting to derive an interpretation for 
the existence of the PHP on the basis of more complex models is well beyond the scope of the 
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present work, and therefore we do not discuss further here the above mentioned hypotheses. 
As detailed in Appendix [HI it can be shown that the specific entropy is given by the following 
formula: 

s = 0.51nlO-logM/L + 2.51nlO-logre+ 

(8) 

- 0.2 InlO </i>c +$(n) + const, 
where ^{n) is an dimensionless function given by: 



P and p being the dimensionless 3D pressure and density profiles of the models respectively. 
Since $(n) is a very complex function of the Sersic index, there is no a priori reason for 
which a systematic variation of s along the galaxy sequence (such as the s-M relation found 
by MLCOO) should imply the existence of a plane in the space of logrg, </i>e and logn. 
From Eq. IHl we see that this is the case only and only if the function is a linear function 
of logn. In Fig. El we plot the function for different values of logra. The values of $(?t,) 
were computed numerically, by using the formulae in Appendix O By looking at the figure, 
it is remarkable to note that there exists almost an exact linear relation between $(n) and 
logn. A linear best-fit gives $(n) ~ 0.988 ■ log(n) + 1.57, with residuals smaller than 0.5% in 
absolute value. We point out that the existence of a linear relation between $(n) and logn 
holds in a wide range of n values, from n ~ 1 to n ~ 10, and it is not a consequence of 
the limited range considered for the Sersic index. We conclude, therefore, that the flatness 
of the PHP of ETGs is strictly connected to the physical origin of this relation, the specific 
entropy of ETGs being a linear combination of the variables logrp, </i>e and logn. In the 
future, it will be very interesting to derive Eq. |Hlby considering more complex galaxy models 
(e.g. by including both dark and luminous matter) and to analyze the implications of such 
models for the existence of a fiat surface in the space of structural parameters. 

7.3 The waveband dependence of the PHP 

Since the variation of </i>e between different wavebands is proportional to the galaxy 
colour (see EqlH} and therefore to galaxy luminosity, through the CM relation, the PHP 
slopes of MS 1008 in the optical and NIR wavebands can be used to constrain how the 
ratios ^ and vary along the galaxy sequence. In Appendix^ we show that, by using 
the i?-band PHP equation and the R-K CM relation, it is possible to derive an equation 




(9) 



+ 2 



?ir(2n) 



■/pin (p3/2^-5/2) ds, 
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Figure 10. Relation between the function <l?(n) and the logarithm of the Sersic index (see the text). The grey curve in the 
upper panel marks the best-fitting line of <I>(7i) vs. logn. The lower panel shows the residuals of the linear fit (solid line). 

which is analogous to that of the K-hand PHP, except for a term which is proportional to 
log It turns out, therefore, that the waveband variation of the PHP does not depend 
significantly on the ratio of effective radii, but it mainly informs on how the ratio of Sersic 
indices vary with other structural parameters of galaxies. 

As detailed in Appendix|Dl assuming that the ratio of shape parameters is only a function 
of the galaxy luminosity, i.e. 

(10) 

where K is the total luminosity in the K band, we obtain two independent constraints on 
uji from the R- and K- band values of a and b respectively. These constraints are 



log = UJiK + UJ2, 

riK 



for the coefficient a, and 



(feK-feR)(i+5feRr) 

aR(l-5feR) ' 



(11) 



(12) 



for the coefficient b, where F is the slope of the R — K CM relation^ of MS 1008. Using the 
values of a and b from Tabled] and Table and taking into account the correlation of uncer- 
tainties on PHP coefficients, we obtain ui = 0.035 ± 0.035 (~ 8 ± 8%) and uji = -0.2 ± 0.2 
(—46 ± 46%) from the first and second equations respectively. We note (i) that the two 

^ For consistency with the procedure of Appendix|D] we calculated F by constructing the CM relation with total galaxy colours, 
obtained from the difference of Kron magnitudes in the R and K bands. This gives F = —0.08 ± 0.02. 
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values of uoi are consistent with each other, showing that the ansatz in Eq. ^Jis compati- 
ble with the values of the PHP slopes, and (ii) that the first equation sets much stronger 
constrains on the value of cji, implying that the variation of log — on galaxy luminosity is 
mainly driven by the logn coefficient of the PHP. Since the value of uJi is fully consistent 
with zero, we conclude that the ratio of optical-NIR shape parameters does not change or 
can have only a mild variation with ga laxy luminosity. This result can be compared with 
the findings of iLa Barbera et al.l (|2003l hereafter LBM03), who showed that the waveband 
variation of the Kormendy relation slope constrains how the ratio of effective radii between 
different wavebands vary with galaxy luminosity. In other terms, the waveband dependences 
of the KR slope and of the PHP logn slope constrain the variation with luminosity of the 
ratios of galaxy radii and shape parameters respectively. This is an interesting result, since 
these ratios fully characterize the radial colour profile in galaxies. Since LBM03 found that 
does not change significantly with luminosity, the above result indicates that all the 
properties of the colour profile in galaxies (that is its shape and gradient) do not change with 
luminosity. The existence of a correl ation between the col our profile and galaxy luminosity 

( 199Cl[ l found n o correlation among the 



is a still debated issue. For example. 



Peletier et al 



Tamura fc Ohta 



(20031) 



internal colour gradients of field ETGs and their luminosities, while 
found that for very bright ETGs in a nearby cluster such correlation could exist. A steep- 
ening of colour gradients with galaxy luminosity is a natural expectation of the monolithic 



collapse model ()L arson 



19741 ) of galaxy formation, since galactic winds blow earlier in less 



massive galaxies, preventing gas dissipation to carry heavy metals in the center, and pro- 
ducing, therefore, a less steep gradient in these systems. On the other hand, in a hierarchical 
scenario of galaxy for mation, merging would dilute stellar population gradients in larger 



galaxies (jWhite 



1980I ). sweeping out or reversing the colour gradient-luminosity relation. 



The implications deriving from a multiwavelength analysis of the KR and the PHP should 
be carefully taken into account by any model aimed to explain the processes underlying the 
formation and evolution of ETGs. 



7.4 The redshift dependence of the PHP 

Since the PHP is a correlation among global parameters of ETGs which are strictly related 
both to their luminosity density (</i>e) and internal structure {r^ and n), its slopes indicate 
how these properties vary along the galaxy sequence. The fact that the PHP slopes do not 
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change significantly up to ^ ~ 0.3 implies that in this redshift range (i) the luminosity 
evolution of (bright) ETGs is not significantly different along the galaxy sequence and (ii) 
correlations among structural properties and galaxy mass do not change significantly with 
z. Point (i) follows directly from the PHP equation (Eq. P). If the evolution with 2; of </i>e 
would depend on log re and/or n, we would expect to find a change of a and b with redshift 
(see Eq. P). The fact that the luminosity evolution of (bright) ETGs does not change along 



the galaxy se quence is consistent with results o : 



Barger et al 



Ziegler et al 



Korniendy relation studies a t intermediate 



199fl : 



La Barbera et al, 



2003), which found 



redshifts (e.g. 

that the luminosity evolution of ETGs is almost independent of galaxy size. Point (ii) can 
be analyzed by considering Eq. |H1 If the PHP origins from a correlation between specific 
entropy and galaxy mass, the slopes of the PHP are fully characterized by the slopes of 
the M/L vs. M and s vs. M relations. Since FP studies seem to indicate that the M/L 



vs. M relation does not change significantly at 2; < 0.3 (jKelson et al 



20001 ). the result of 



Section in21 implies that the slope of the s vs. M relation does not change significantly with 
z. This is consistent with the idea that merging is the physical process which builds up the 



s-M relati on (MLCOO) 



Kauffmann 



and that bright ETGs in clusters are mostly assembled at redshift 



19951). 



z>l (e.f 

As shown in Section 16.21 the zero-point of the PHP can be applied to constrain the 
mean luminosity evolution of ETGs with redshift. This use of the PHP is based on the 
fact that its slopes do not evolve with redshift and that the variations with z of the mean 
values of both log re and logra are known exactly (see Eq. Assuming A(logre) = and 



A(logn) = 0, we find that the evolution of the zero-point of the PHP from z ~ 0.3 to 2 ~ 
is consistent with the cosmological dimming of mean surface brightness and the passive 
fading of an old stellar population with a high formation redshift, Z{ > 1-2. The hypothesis 
A(logre) = is well motivated by the findings of different studies tha t the distribution 



of radii of ETGs does not change significantly at intermediate redshifts ()La Barbera et al 



mm 



Shen et al 



20031 ) . On the other hand, assuming A(logri) = is a more tricky point. 



since a morphq 



(Dressier et al. 



■Ogv-d ensitv relation is known to exist and this relation evolves with redshift 



19971 ). To account for this effect, we calculated the mean values of logra for 



the sample of MS 1008 and that of GRA02. The difference between the two values turned 
out to be A(log?2) = —0.045 ± 0.05^, in agreement with the hypothesis A(logn) = 0. 



This value was obtained by considering only the galaxies of the GRA02 sample within the same magnitude range of the 
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Due to heavy request for measuring velocity dispersions of distant galaxies, it is clearly 
of great interest to have shown that the PHP is a valuable tool for measuring the luminosity 
evolution of ETGs. In the future, it will be very interesting to derive the PHP for samples of 
galaxies at higher redshift, taking advantage (a) of the relatively small dispersion of the PHP 
with respect to other purely photometric correlations, and (b) of the fact that the PHP slopes 
seem not to depend on the waveband, allowing, therefore, a straightforward comparison of 
the PHP among different redshifts (different restframe bands) to be performed. 

8 SUMMARY 

Using a unique dataset of optical and NIR data, we have done a detailed study of the 
Photometric Plane (PHP) relation in the i?, / and K bands for a large sample of ETGs 
in the rich galaxy cluster MS 1008 at intermediate redshift {z ~ 0.3). We have derived the 
PHP by accounting for selection effects and other issues related to the fit of this relation. 
Our main findings can be summarized as follows. 

(i) The PHP is indeed a plane in the space of the quantities logrg, <fi>e and logn, 
i.e., there is no significant departure from a purely flat relation (Section I5.1|l . We have 
shown that this result is consistent with the fact that the PHP origins from a systematic 
variation of specific entropy of ETGs along the galaxy sequence. For isotropic, non rotating, 
one component models, the specific entropy turns out, in fact, to be almost an exact linear 
combination of logrc, </i>e and logn. 

(ii) We have found no waveband dependence of the PHP slopes f Section 16.11) . Both the 
PHP coefficients and its intrinsic scatter are fully consistent among the optical and NIR 
wavebands. 

(iii) The scatter around the PHP is about ~ 32 %, which is larger than the average 
scatter of the Fundamental Plane, but smaller than that of the Kormendy relation and the 
photometric or the 'Mg2 Fundamental Plane' (Section 7.1). 

By comparing the PHP coefficients at z = 0.3 with those at z ~ for the sample of 
GRA02, we have found that the slopes of the PHP do not show a significant variation 
in this redshift range. This fact has important consequences on our understanding of the 
evolution of bright ETGs, as discussed in Section 7.4. Namely: 

galaxies in MS 1008 (i.e. ^ m * +3). The same result, however, is obtained by not applying any magnitude cut. In this case we 
obtain A (log n) = -0.04 ± 0.06 
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(iv) the luminosity evolution does not change significantly along the galaxy sequence; 

(v) the slope of the relation between specific entropy and galaxy mass does not change 
since z ~ 0.3. 

Both these findings are in agreement with the fact that bright ETGs are already assembled 
at high redshift. Finally, we have demonstrated that the PHP is a valuable tool for measuring 
the luminosity evolution of ETGs. We have found that 

(vi) the mean luminosity of ETGs is consistent with the passive fading of an old stellar 
population, with formation redshift > 1. 

Due to the heavy request of telescope time for measuring velocity dispersions, construct- 
ing the FP relation for large samples of galaxies at 2; > 0.5 becomes impracticable. The 
present results show that the PHP could be a very interesting alternative tool for the study 
of ETGs at high redshifts. 
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APPENDIX A: PSF MODELING 

The stars in the R, I, K{S) and K{P) images were fitted by the following formula: 
PSF(r) = X:Ak(0)-Mk(r,^), (Al) 



k=i 



where r is the distance to th e star center, 6 is the polar angle, M oc [1 + (r/rc)^] ^ is 



the Moffat law ( Moffat 



1969| ). p is the number of Moffat functions, and A is an angular 



modulation function. PSF asymmetries were described by adopting the formula: 



30 F. La Barbera et al. 

A{e) = ^ [a„, cos(m ■ 9) + sm(m ■ 9)] . (A2) 

m=l 

The values of p and q were chosen interactively for each image in order to obtain an accurate 
modeling of stellar isophotes. This was achieved by using two or three Moffat functions, 
depending on the image, with g ^ 3. To account for PSF variations across the frames, each 
galaxy was fitted by using a local PSF model, obtained from the nearest star in the field. The 
PSF modeling in the R band image is illustrated in Fig. lAH where we show residuals of star 
fitting as a function of both the distance to the star center and the polar angle. Residuals 
are plotted for two fitting cases, with and without the angular modulation function. The 
figure shows that the adopted models give an accurate description of the PSF, allowing both 
the radial and the angular behaviour of the PSF to be reproduced. We note that the use 
of asymmetric terms in the fit allows a better centering of the PSF model to be obtained. 
Hence the improvement in the radial trend of PSF residuals which is observed in Fig. lAll 
Similar results were also obtained in the analysis of the images of MS 1008 in the other 
bands. 

APPENDIX B: UNCERTAINTIES ON r^, <fi>^ AND n 

Uncertainties on structural parameters were estimated by applying the 2D fitting method to 
simulated galaxy images. For each galaxy in the R, I, K{S) and K{I) bands, we used the 
best-fitting values of rg, < /i >e and n to construct a set of seeing-convolved Sersic models to 
which both photon and read-out noise were added. In order to account for uncertainties due 
to PSF modeling and to image masking, each simulation was obtained by using a different 
PSF, according to the uncertainties on the PSF fitting parameters, while the 2D fitting 
was performed by using the same PSF and the same mask of the real galaxies. For each 
galaxy, uncertainties on structural parameters were estimated from the covariance matrix of 
2D best-fitting parameters for the corresponding set of simulations. The mean uncertainties 
on logTg, </i>e and logn amount to ~ 0.03 dex, ~ 0.18 mag/arcsec^ and ~ 0.05 dex 
respectively, in R; to ~ 0.05 dex, ~ 0.26 mag/arcsec^ and ~ 0.09 dex in /; to ~ 0.06 dex, 
~ 0.27 mag/arcsec^ and ~ 0.06 dex in K{I); and to ~ 0.07 dex, ~ 0.34 mag/arcsec^ and 
~ 0.08 dex in K{S). 

In Fig. IBll we compare the structural parameters of the = 1 1 2 galaxies in common 
between the R and I bands. We note that the mean differences of the distributions are fully 
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Figure Al. PSF fitting to four stars in tlie FORSl ii-band image. Eacli plot shows tlie image of the star (upper-left panel), 
the residual map (upper-right) , and mean fitting residuals in units of the background standard deviation as a function of the 
radial and angular coordinates (lower-left and lower-right panels, respectively). Large gray symbols and small black circles 
denote the residuals obtained by respectively excluding and including the angular modulation function in the PSF fit (see text 
for details). The spatial scale of each image is the same of that shown in the upper-left plot. 



consistent with zero, in agreement with the fact that R and / bands at 2; ~ 0.3 sample 
a very similar spectral region and ETGs have small optical-optical colour gradients. For 
this reason, the widths of the distributions in Fig. IBll provide a rough estimate of the mean 
uncertainty on logrg, </i>e and logn. These values are fully consistent with those obtained 
by adding in quadrature the above mean uncertainties on structural parameters. In Fig. IB2I 
we compare the structural parameters of the = 38 galaxies in common between the -^"(5") 
and the K{I) samples. Again, the mean differences of the distributions are fully consistent 
with zero, showing that no significant systematic effects are present. We also note that the 
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m = 0.003 
5(m) = 0.006 
std = 0.06 



' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' ' ' ' I ' 
m = -0.013 
5(m) = 0.024 
std = 0.24 



m = -0.007 
5(m)=0.007 
std = 0.07 
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A(log r^) A(</.t>^) A(log n) 

Figure Bl. Comparison of structural parameters for the N = 112 galaxies in common between the R- and 7-band samples. The 
quantities m, <5(m) and std are the mean value, the uncertainty on the mean, and the standard deviation for each distribution. 
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Figure B2. The same of Fig.l^for the = 38 galaxies in common between the K(S) and K{I) samples. 



standard deviations shown in Fig. IB2I agree with what expected on the basis of the mean 
errors on structural parameters (see above). 



APPENDIX C: SPECIFIC ENTROPY OF ETGS 

As shown by MLCOO, the specific entropy of a galaxy, s, is given by the following formula: 

s = 1 1 (in p3/2^-V2^ p dV, (CI) 

where p and P are the 3D density and pressure fields, respectively, while Mt is the total 
mass. Following LGM99 and MLCOO, we consider spherical, non-rotating galaxy models, 
with negligible radial gradients of the mass-to-light ratio. For such models the function p(r) 
is obtained directly from the 2D Sersic law by solving the Abel integral equation, while P{r) 
is calculated from the equation of hydrostatic equilibrium, ^ = —GMp/r'^, where the mass 
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profile M(r) is obtained from a direct integration of p(r). The functions p(r), M(r) and 
P(r) can be written as products of dimensional factors for dimensionless functions: 

p(r) = M/L ■ lo/Re ■ p{s) (C2) 
M(r) = M/L -lo-Rl- M{s) (C3) 
P(r) = G{M/L ■ lof ■ P{s) (C4) 

where s = r/R^ is the dimensionless radius and Jq is the 2D central surface brightness. 
By substituting the previous formulae in Eq. ICll and by using the relation between Jq and 
</>e (Eq. El), we obtain Eqs. IHlandU 

APPENDIX D: OPTICAL-NIR PHP 

Considering the optical-NIR colour-magnitude relation, R — K = T- K + ^, where R — K 
and K are total galaxy colours and magnitudes, and using the definition of total magnitude 
(Eq. ini), we can write the R-band PHP relation, logi?e,R = ■ logriR + 6r- < /i >e,R +cr, 
where the subscript R denotes the waveband, as follows: 

l0gt; + l0gi?e,K = 

Or log ^ + Or log nK + 6r (< p >c,R - < p >c,k) + 

+6r <p>c,K +Cr = (Dl) 
= Or log ^ + Or log + 6r [P < p >c,K "SP log i?e,K + 

-2.5P log(27r) + A + 5 log fj] + 6r < /i >e,K +cr. 
From this equation, we obtain: 

logPe,K = (1 + 56rP)-^ • {aRlog ^ + (56r - 1) log 1^+ ^^^^ 
aRlognK + &R(l + P) <Ai>c,K + [cR + feRA-2.5Plog(27r)]}. 

We note that since (56r — 1) ~ (see TableE}, the term log in Eq. lD2l can be neglected. 
Moreover, since (1 + 56rP)~^ ■ (1 + P) ~ 1, we obtain from the previous equation: 

logi?e,K ^ (1 + 56rP)-^- aRlog ^ + (1 + 56rP)-^- 

■Or log nK + &R < /i >c,K +const. 
We note that this equation is similar to that of the PHP in the K band, except for a 
term which is proportional to the ratio of shape parameters between the R and K bands. 
Substituting the ansatz of Eq. Elin the previous equation, we obtain a relation identical to 
that of NIR PHP. Comparing the slopes of this relation with the values of a and h in the K 
band, we obtain Eqs. [TT] and [T21 respectively. 
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1 n 




— 12 


41 




-0.53±0.08 


19.48±0.46 


0.70±0.07 


-0.52±0.18 


18.84±1.05 


0.67±0.10 


-0.65 ±0.02 


15.51 ±0.08 


0.49±0.04 


in 


1 n 

lU 


OQ fifi 


1 

— iz 


41 




-0.49±0.02 


19.92±0.13 


0.46±0.03 


-0.49±0.01 


19.26±0.05 


0.46±0.05 


-0.54±0.01 


16.25±0.05 


0.32±0.03 


1 n 


1 n 

lU 


OO AC 
ZZ. Do 


1 o 
— iz 


/1 1 
^1 


UU.O 


0.00±0.14 


22.00±0.60 


0.79±0.10 


0.01 ±0.03 


21.35±0.23 


0.81±0.13 








in 


1 n 

lU 


OR '7A 
Zo. / "+ 


— 12 


41 


UU. 1 


-0.57±0.05 


18.87±0.37 


0.58±0.06 


-0.57±0.16 


18.13±1.20 


0.72±0.27 


-0.48±0.01 


15.98±0.04 


0.33 ±0.02 


lU 


1 ri 

lU 




1 o 

— IZ 


/in 


zu.u 


-0.31 ±0.06 


19.52±0.42 


0.78±0.14 


-0.26±0.17 


19.32±0.78 


0.60±0.06 


-0.34±0.03 


16.14±0.16 


0.74±0.07 


1 n 


1 n 

lU 


OR 3*^ 


— 12 


zin 


44 8 


-0.49 ±0.04 


20.88±0.25 


0.30±0.04 














1 n 


1 n 

lU 


•^O RQ 


1 o 

— iz 


/in 


Zl.Z 


-0.11 ±0.06 


20.76±0.24 


0.67±0.07 


-0.13±0.18 


19.99±0.78 


0.63±0.15 


-0.23 ±0.02 


16.62±0.08 


0.64±0.02 


1 u 


1 u 




1 o 

— 1 z 


/in 


JZ. / 


-0.16 + 0.03 


21.60±0.16 


0.44 + 0.03 


-0.13 + 0.03 


21.03 + 0.20 


0.49 + 0.03 


-0.16 + 0.02 


18.13±0.09 


0.44+0.03 








1 o 
— IZ 


/in 


z.y).z. 


0.31±0.1() 


22.17±0.47 


0.88 + 0.06 


0.23±0.06 


21.13±0.20 


0.88±0.16 


-0.05±0.09 


16.91±0.35 


0.87±0.06 


1 n 


1 n 


1A 7^ 


— 12 


40 


27 1 


-0.58±0.02 


20.26 ±0.1 6 


0.35 ±0.09 


-0.54±0.38 


19.81±1.70 


0.61 ±0.12 








1 n 


1 n 

lU 


30 1 7 


— 12 


An 




-0.38±0.02 


20.66±0.08 


0.59±0.04 


-0.38±0.06 


20.05 ±0.51 


0.59±0.12 


-0.47±0.03 


16.97±0.13 


0.56±0.08 


1 n 


1 n 

lU 


OQ fl7 


1 o 

— iz 


/in 


1 ^ Q 

ID. J 


-0.49±0.04 


19.98±0.24 


0.45±0.03 


-0.50±0.05 


19.25±0.45 


0.48±0.12 


-0.49±0.02 


16.50±0.07 


0.47±0.02 


1 n 

lU 


1 n 


OQ 


1 o 

— IZ 


Afi 




-0.23 ±0.04 


20. 39 ±0.21 


0.56±0.02 


-0.22±0.04 


19.76±0.22 


O.55±0.03 


-0.28±0.01 


16.57±0.06 


0.53±0.02 




lU 


A ^ OA 


1 o 
— IZ 


/in 




-0.32±0.06 


20.47±0.35 


0.58±0.04 


-0.13±0.06 


20.38±0.43 


0.96±0.10 








1 n 


1 n 

lU 


3il 3il 


— 12 


An 


no 4 


-0.22±0.03 


20.82±0.19 


0.53±0.03 


-0.21 ±0.04 


20.18±0.23 


0.55±0.03 


-0.21 ±0.02 


17.22±0.10 


0.63±0.03 


1 n 


1 n 

lU 


07 7A 
Z /. /O 


— 12 


An 


1 J. / 


-0.81 ±0.57 


19.31 ±3.34 


0.54±0.20 








-0.98±0.04 


14.81 ±0.19 


0.28±0.07 


10 


10 


28.07 


-12 


40 


09.6 


-0.48±0.03 


20.18±0.15 


0.67±0.07 


-0.44±0.07 


19.66±0.58 


0.63±0.20 


-0.48±0.03 


16.83±0.14 


0.51 ±0.05 


10 


10 


33.52 


-12 


40 


04.6 


-0.64 ±0.02 


19.27±0.09 


0.51 ±0.03 


-0.61 ±0.08 


18.67±0.61 


0.52±0.09 


-0.74±0.02 


15.14±0.09 


0.53±0.05 


10 


10 


33.45 


-12 


40 


02.4 


-0.82±0.28 


18.94±1.68 


0.50±0.11 








-1.02±0.05 


14.32±0.29 


0.53±0.12 


10 


10 


37.43 


-12 


39 


56.4 


-0.18±0.04 


22.45±0.16 


0.42±0.06 


0.15±0.37 


22.93 ±1.36 


0.62±0.28 








10 


10 


42.09 


-12 


39 


45.4 


-0.14±0.05 


21.03±0.25 


0.69±0.06 


-0.17±0.06 


20.30±0.22 


0.57±0.13 


-0.22±0.04 


17.21 ±0.19 


0.58±0.06 


10 


10 


34.31 


-12 


39 


52.0 


-0.52±0.02 


19.67±0.13 


0.49±0.03 


-0.50±0.06 


19.03 ±0.47 


0.50±0.06 


-0.61 ±0.16 


15.47±0.86 


0.71 ±0.06 


10 


10 


46.01 


-12 


39 


41.1 


0.09 ±0.05 


21.78±0.25 


0.51 ±0.03 


0.12±0.24 


21.14±0.84 


0.57±0.28 








10 


10 


29.84 


-12 


39 


45.3 


-0.15±0.02 


20.88±0.14 


0.35±0.03 


-0.14±0.03 


20.24±0.17 


0.40±0.03 


-0.14±0.01 


17.35 ±0.07 


0.48±0.02 


10 


10 


29.45 


-12 


39 


50.0 


-0.68±0.07 


18.62±0.42 


0.30±0.02 








-0.64±0.01 


15.53 ±0.05 


0.45±0.02 


10 


10 


29.09 


-12 


39 


35.6 


-0.19±0.05 


19.96±0.26 


0.73±0.04 


-0.13±0.06 


19.64±0.39 


0.74±0.13 


-0.32±0.23 


15.81 ±1.05 


0.77±0.05 


10 


10 


31.08 


-12 


39 


43.3 


-0.12±0.03 


21.32±0.15 


0.41 ±0.03 


-0.26±0.04 


19.97±0.26 


0.52±0.03 


-0.32±0.01 


16.82±0.06 


0.60±0.02 


10 


10 


27.01 


-12 


36 


38.5 


-0.62±0.04 


19.22±0.24 


0.61 ±0.04 


-0.54±0.15 


18.90±0.83 


0.50±0.21 








10 


10 


26.59 


-12 


36 


38.1 


-0.48±0.01 


19.92±0.10 


0.37±0.03 


-0.46±0.01 


19.33 ±0.09 


0.32±0.02 








10 


10 


33.03 


-12 


37 


27.6 


-0.67±0.03 


18.87±0.17 


0.58±0.03 


-0.64±0.14 


18.49 ±1.06 


0.48±0.22 








10 


10 


32.53 


-12 


37 


04.6 


0.03±0.08 


21.95±0.36 


0.75±0.05 


0.03±0.03 


21.28±0.22 


0.72±0.13 








10 


10 


45.21 


-12 


39 


13.1 


-0.27±0.11 


20.80±0.53 


0.82±0.08 


-0.29±0.09 


20.02 ±0.63 


0.83±0.12 


-0.53 ±0.39 


16.17±1.76 


0.83±0.11 


10 


10 


26.93 


-12 


37 


07.8 


-0.37±0.03 


20.21 ±0.20 


0.52±0.03 


-0.38±0.04 


19.44±0.33 


0.58±0.08 
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RAj2000 


DECj2ooo 


logre,R 


<M>e,R 


lognR 


log r^.i 


<A'>e,I 


logni 


logre.K 


<A'>e.K 


log njT 














(") 


(magl"^) 




(") 


(magl"^) 




(") 






10 


10 


32.33 


-12 


37 


19.4 


0.18±0.06 


21.75±0.27 


0.59±0.03 


0.16±0.02 


21.01±0.12 


0.57±0.03 


0.04±0.02 


17.70±0.10 


0.49±0.03 


10 


10 


32.87 


-12 


36 


59.5 


-0.26±0.04 


19.93±0.15 


0.69±0.06 


-0.27±0.07 


19.24±0.33 


0.60±0.11 








10 


10 


36.43 


-12 


37 


15.9 


0.42±0.08 


22.54±0.33 


0.64±0.02 


0.61±0.13 


22.55±0.51 


0.71 ±0.05 








10 


10 


22.29 


-12 


37 


26.7 


-0.41 ±0.03 


19.77±0.20 


0.47 ±0.02 


-0.36±0.04 


19.32±0.38 


0.47 ±0.1 9 








10 


10 


33.96 


-12 


39 


14.3 


-0.15±0.08 


21.25±0.35 


0.87±0.O8 


-0.20±0.17 


20.39±0.82 


0.85±0.14 


-0.13±0.15 


17.99±0.72 


0.90±0.08 


10 


10 


27.94 


-12 


37 


37.7 


0.22±0.10 


22.93±0.40 


0.59±0.06 


0.26±0.06 


22.44±0.26 


0.61 ±0.03 


0.08±0.17 


19.21 ±0.62 


0.42±0.14 


10 


10 


31.24 


-12 


37 


48.5 


-0.04 + 0.05 


20.36±0.25 


0.58±0.02 


-0.03 ±0.04 


19.73 + 0.21 


0.57±0.03 








10 


10 


33.53 


-12 


37 


46.7 


0.05 ±0.04 


22.45±0.18 


0.35 ±0.04 


0.06±0.07 


22.02 + 0.20 


0.36±0.14 








10 


10 


28.92 


-12 


37 


45.3 


-0.52 ±0.06 


20.06±0.35 


0.70 ±0.07 


-0.52 ±0.11 


19.41+0.85 


0.67±0.17 


-0.45±0.04 


17.01 ±0.19 


0.44±0.09 


10 


10 


43.98 


-12 


39 


06.9 


-0.51 ±0.05 


21.08±0.20 


0.57±0.09 


-0.50±0.05 


20.49 +0.40 


0.52±0.06 








10 


10 


34.11 


-12 


37 


46.0 


-0.64±0.04 


20.31 ±0.22 


0.33±0.09 


-0.60±0.03 


19.86+0.13 


0.31 ±0.06 








10 


10 


37.78 


-12 


39 


10.8 


-0.31±0.05 


20.94±0.24 


0.71 ±0.04 


-0.35±0.09 


20.17±0.65 


0.66±0.13 


-0.48±0.10 


16.81±0.44 


0.62±0.07 


10 


10 


29.08 


-12 


39 


14.9 


-0.29±0.10 


21.30±0.55 


0.81±0.12 


-0.44±0.17 


19.84±1.28 


0.89±0.31 


-0.55 ±0.08 


16.56±0.45 


0.85±0.09 


10 


10 


42.21 


-12 


38 


55.5 


0.63±0.11 


22.99±0.43 


0.74±0.02 


0.70±0.19 


22.55±0.75 


0.86±0.13 


0.50±0.13 


18.80±0.53 


0.90±0.06 


10 


10 


36.62 


-12 


37 


51.2 


-0.47±0.02 


19.98±0.12 


0.31±0.02 


-0.44±0.04 


19.41 ±0.24 


0.32±0.10 








10 


10 


33.01 


-12 


38 


53.1 


-0.70±0.02 


19.14±0.13 


0.30±0.04 


-0.68±0.14 


18.59±0.85 


0.38±0.19 


-0.67±0.06 


15.99±0.34 


0.42±0.12 


10 


10 


32.97 


-12 


38 


55.8 


-0.12 + 0.04 


21.57±0.19 


0.42 + 0.05 


-0.16 + 0.02 


20.82 + 0.13 


0.34±0.01 


-0.04±0.05 


18.53±0.19 


0.49±0.04 


10 


10 


33.10 


-12 


37 


57.8 


-0.59±0.02 


19.62±0.11 


0.39+0.03 


-0.58 + 0.04 


18.99+0.35 


0.37±0.10 








10 


10 


30.80 


-12 


37 


44.8 


-0.71 ±0.02 


19.17±0.14 


0.37±0.02 


-0.67+0.05 


18.69+0.22 


0.37±0.06 








10 


10 


25.86 


-12 


38 


10.4 


-0.05±0.08 


21.21 ±0.34 


0.66±0.07 


-0.03±0.05 


20.64±0.27 


0.67±0.02 


-0.17±0.02 


17.37±0.08 


0.54±0.03 


10 


10 


28.59 


-12 


38 


18.9 


-0.17±0.04 


20.96±0.21 


0.64±0.05 














10 


10 


27.87 


-12 


38 


14.3 


-0.74 ±0.03 


19.87±0.20 


0.35 ±0.08 


-0.73±0.05 


19.23±0.42 


0.35±0.07 








10 


10 


30.02 


-12 


38 


18.5 


-0.21 ±0.04 


20.64±0.19 


0.63 ±0.03 


-0.23±0.05 


19.89±0.34 


0.61 ±0.08 


-0.27±0.03 


16.96±0.13 


0.49±0.06 


10 


10 


27.45 


-12 


38 


22.9 


-0.73±0.03 


19.91±0.16 


0.39±0.08 


-0.71 ±0.08 


19.33±0.61 


0.44±0.08 








10 


10 


33.31 


-12 


38 


35.3 


-0.04 + 0.08 


21.11±0.36 


0.69±0.05 


-0.08±0.04 


20.32±0.29 


0.67±0.08 


-O.12±0.03 


17.38±0.15 


0.63 + 0.06 


10 


10 


33.12 


-12 


38 


35.8 


-0.30 + 0.04 


19.92±0.21 


0.62±0.02 


-0.28±0.05 


19.32±0.47 


0.64±0.17 


-0.34±0.03 


16.14±0.15 


0.58 + 0.06 


10 


10 


29.98 


-12 


38 


36.0 


-0.48 ±0.02 


19.99±0.09 


0.59 ±0.05 


-0.50±0.07 


19.26±0.51 


0.53±0.07 


-0.48±0.03 


16.61±0.16 


0.48 + 0.07 


10 


10 


38.82 


-12 


38 


22.4 


-0.24±0.05 


20.92±0.22 


0.73 ±0.06 








-0.44±0.18 


16.41 ±0.90 


0.75+0.16 


10 


10 


27.35 


-12 


38 


55.5 


-0.45±0.02 


20.65±0.10 


0.59±0.04 


-0.45±0.08 


20.03 ±0.58 


0.50±0.25 


-0.53 ±0.03 


16.96±0.12 


0.41+0.05 


10 


10 


31.16 


-12 


38 


41.0 


-0.32±0.02 


20.00±0.13 


0.59±0.02 


-0.39±0.04 


19.07±0.40 


0.55±0.16 


-0.41 ±0.01 


16.12±0.07 


0.52 + 0.04 


10 


10 


31.27 


-12 


38 


46.7 


-0.80±0.09 


18.99±0.60 


0.39±0.11 


-0.79±0.05 


18.35±0.49 


0.38±0.11 


-0.85 ±0.21 


15.21 ±1.05 


0.49+0.11 


10 


10 


30.49 


-12 


38 


36.5 


-0.22±0.13 


22.86±0.47 


0.41±0.19 














10 


10 


46.13 


-12 


38 


32.0 


-0.20±0.03 


21.68±0.16 


0.33±0.03 


-0.18±0.02 


21.14±0.10 


0.34±0.08 








10 


10 


26.69 


-12 


38 


40.9 


-0.65±0.02 


20.06±0.17 


0.31 ±0.07 


-0.65±0.05 


19.41 ±0.43 


0.31 ±0.06 








10 


10 


32.23 


-12 


38 


42.1 


-0.71 ±0.04 


20.44±0.24 


0.36±0.15 














10 


10 


28.82 


-12 


38 


44.6 


-0.61 ±0.05 


20.46±0.27 


0.52±0.09 


-0.61 ±0.09 


19.85 ±0.70 


0.49±0.14 








10 


10 


27.78 


-12 


38 


47.2 


-0.71 ±0.04 


20.06±0.26 


0.31±0.22 


-0.68±0.10 


19.53 ±0.62 


0.34±0.12 








10 


10 


28.20 


-12 


38 


52.0 


-0.64±0.04 


19.68±0.25 


0.32±0.05 


-0.64±0.35 


19.03 ±2.07 


0.31 ±0.40 


-0.61 ±0.01 


16.35±0.05 


0.31+0.04 


10 


10 


34.74 


-12 


39 


01.2 


-0.17 ±0.04 


21.50±0.20 


0.63 ±0.05 


-0.23±0.03 


20.61 ±0.22 


0.58±0.05 


-0.22±0.10 


17.74±0.39 


0.64+0.06 


10 


10 


34.16 


-12 


38 


56.3 


-0.68 + 0.03 


19.57±0.17 


0.53±0.06 


-0.69 ±0.04 


18.84±0.39 


0.56±0.06 


-0.73±0.10 


15.96±0.56 


0.65 + 0.08 


10 


10 


36.90 


-12 


39 


00.5 


-0.60+0.02 


20.32±0.14 


0.44±0.04 


-0.57±0.07 


19.77±0.52 


0.44±0.11 


-0.55 ±0.08 


16.97±0.32 


0.50+0.09 


10 


10 


32.79 


-12 


39 


09.5 


-0.08±0.05 


20.78±0.24 


0.71 ±0.06 


-0.13±0.04 


19.90±0.25 


0.65±0.03 


-0.25±0.02 


16.58±0.10 


0.66+0.03 


10 


10 


27.67 


-12 


39 


00.8 


-0.44±0.03 


20.50±0.18 


0.49±0.04 


-0.45±0.02 


19.84±0.09 


0.46±0.08 


-0.52±0.02 


16.88±0.09 


0.39+0.04 


10 


10 


35.12 


-12 


39 


07.1 


-0.34±0.04 


20.96±0.18 


0.67 ±0.06 


-0.35±0.08 


20.25 ±0.60 


0.63±0.14 


-0.45 ±0.03 


17.03±0.12 


0.60+0.04 


10 


10 


31.17 


-12 


39 


10.0 


-0.47±0.06 


20.44±0.27 


0.73 ±0.07 


-0.36±0.12 


20.21 ±0.89 


0.87±0.24 


-0.37±0.09 


17.29±0.40 


0.92+0.06 


10 


10 


22.92 


-12 


39 


20.6 


-0.32±0.08 


21.48±0.37 


0.57±0.09 


-0.29±0.05 


21.02±0.38 


0.58±0.10 








10 


10 


34.81 


-12 


39 


20.8 


-0.33±0.03 


22.10±0.13 


0.36±0.06 


-0.30±0.02 


21.56±0.11 


0.38±0.06 








10 


10 


18.71 


-12 


39 


33.9 


-0.45 ±0.03 


19.76±0.17 


0.39±0.02 


-0.42±0.06 


19.23±0.43 


0.38±0.20 








10 


10 


32.98 


-12 


39 


26.7 


-0.29±0.01 


21.29±0.09 


0.30±0.04 


-0.25±0.02 


20.76±0.18 


0.36±0.16 


-0.17±0.04 


18.22±0.15 


0.54+0.03 


10 


10 


31.17 


-12 


39 


31.4 


0.06±0.06 


21.81±0.32 


0.77±0.06 
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